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Dynamics of Lattice Triangulations on Thin Rectangles 

Pietro Caputo, Fabio Martinelli* Alistair Sinclairl Alexandre Stauffer^ 


Abstract 

We consider random lattice triangulations of n x fc rectangular regions with weight where 
A > 0 is a parameter and \a\ denotes the total edge length of the triangulation. When A S (0,1) 
and k is fixed, we prove a tight upper bound of order v? for the mixing time of the edge-flip 
Glauber dynamics. Combined with the previously known lower bound of order exp(0(n^)) for 
A > 1 |3], this establishes the existence of a dynamical phase transition for thin rectangles with 
critical point at A = 1. 


1 Introduction 


Consider an n x /c lattice rectangle ^ = {0,1,..., n} x {0,1,..., /c} in the plane. A triangulation 
of is defined as a maximal set of non-crossing edges (straight line segments), each of which 
connects two points of and passes through no other point. See Figure 1 for an example. 




Figure 1: Two triangulations of a 5 x 3 rectangle 

Call Q{n,k) the set of all triangulations of All a G Tl{n,k) have the same number of edges 

and the set of midpoints of the edges of a does not depend on a. Thus, we may view a G n(n, k) 
as a collection of variables {ax, x G where 

^n,k • 2’ ^ {d,2)l-)2’'''T 2 T} \ ^n,ki 

is the set of all midpoints. Moreover, any element a G kl{n,k) is unimodular, i.e., each triangle in 
a has area see, e.g., BUM for these standard structural properties. If an edge ax of a is the 
diagonal of a parallelogram, then it is said to be flippable: one can delete this edge and add the 
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opposite diagonal to obtain a new triangulation o' G Q(n,k). In this case o,o' differ by a single 
diagonal flip and are said to be adjacent. The corresponding graph with vertex set and 

edges between adjacent triangulations, called the flip graph, is known to be connected and to have 
interesting structural properties; see PE] and references therein. 

We consider the following model of random triangulations. Fix A G (0, oo) and define a probability 
measure // on k) by 

AkI 

where Z = Yla'en{n k) 1^1 total ii length of the edges in o, i.e., the sum of the 

horizontal and vertical lengths of each edge. The case A = 1 is the uniform distribution, while 
A < 1 (respectively, A > 1) favors triangulations with shorter (respectively, longer) edges. We refer 
to [3] and references therein for background and motivation concerning this choice of weights. 

A natural way to simulate triangulations distributed according to /r is to use the edge-flip Glauber 
dynamics defined as follows. In state o, pick a midpoint x G ^n,k uniformly at random; if the 
edge Ox is flippable to edge cj(, (producing a new triangulation o'), then flip it with probability 

p,{o') + pi{o) + Al‘^^^1 ’ 

else do nothing. Since the flip graph is connected, this defines an irreducible Markov chain 
on n(n, k), and the flip probabilities Q ensure that the chain is reversible with respect to n- Hence 
the dynamics converges to the stationary distribution /r. We analyze convergence to stationarity 
via the standard notion of mixing time, defined by 

Tmix = inf {t G N : max ||p*(fT, •) - fi\\ < 1/4} , (2) 

a&n{n,k) 

where p'{o,-) denotes the distribution after t steps when the initial state is o, and ||z/ — /i|| = 
5 Yla£n{n k) l^(^) “ h{^)\ is the usual total variation distance between two distributions fi, v. 

As discussed in [3|, there is empirical evidence that the value A = 1 represents a critical point 
separating the sub-critical regime A G (0,1), characterized by rapid decay of both equilibrium and 
dynamical correlations, from the super-critical regime A > 1, characterized by the emergence of long- 
range correlations and a dramatic slowdown in the convergence to equilibrium. We substantiated 
this picture by showing that there exist constants C > 0 and Ai G (0,1) such that 

Fmix < Ckn{k -\- n), 

for all k,n € N and for all A < Ai; see P Theorem 5.1]. This estimate is based on a coupling 
argument that requires A to be sufficiently small; in particular, Ai = 1/8 suffices. We conjectured 
in P that the mixing time should satisfy Tmix = 0{kn{k -\- n)) throughout the sub-critical regime 
A G (0,1). However, except for the special case k = 1, establishing even an arbitrary polynomial 
bound on Tmix in the whole region A < 1 has turned out to be very challenging. Regarding the 
super-critical regime, by P Theorem 6.1 and Theorem 6.2] it is known that, for A > 1, one has 
Tmix = exp(H(/c -|- n)) for all k, n, and that = exp(H(n^/A:)) if n > . 

In this paper we establish the conjectured behavior for all A < 1 in the case of “thin” rectangles, 
i.e., the case when k is fixed and n is large. 

Theorem 1.1. For any A G (0,1), /c G N, there exists a constant C = C{\,k) > 0 such that the 
mixing time of the Glauber dynamics for n x k triangulations satisfies Tmix <Cn^ for all n>l. 
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We remark that the above bound is sharp up to the value of the constant C since it is known 
that Tmix > Cokn{k + n) for some positive constant Co for any k,n £ N and any A > 0; see [31 
Proposition 6.3]. However, as a function of k the constant C in Theorem o can be exponentially 
large, and thus the interest of this bound is limited to the case of thin rectangles. 


In the special case k = 1, the above theorem can be obtained by a direct coupling argument; see jH 
Theorem 5.3]. Moreover, it is interesting to observe that in the case k = 1 the set of triangulations 
is in 1-1 correspondence with the set of configurations of a lattice path, and that diagonal flips are 
equivalent to so-called mountain/valley flips in the lattice path representation. Weighted versions 
of lattice path models have been studied extensively in the past (see, e.g., mm), and it is tempting 
to analyze the n x k triangulation model as a multi-path system with k interacting lattice paths. 
While this can be done in principle, it turns out that the interaction between the paths is technically 
very complex. Even the case k = 2 apparently does not allow for significant simplification with this 
representation. 


The proof of Theorem o will rely crucially on some recent developments by one of us m based 
on a Lyapunov function approach to the sub-critical regime A G (0,1). As detailed in subsequent 
sections, the main results of m will be used first to show that after T = 0{n‘^) steps of the chain we 
can reduce the problem to a restricted chain on a “good” set of triangulations, each edge of which 
never exceeds logarithmic length, and then to show that distant regions in our thin rectangles can 
be decoupled with an exponentially small error. This will enable us to set up a recursive scheme 
for functional inequalities related to mixing time such as the logarithmic Sobolev inequality. The 
recursion, based on a bisection approach for the relative entropy functional inspired by the spin 
system analysis of mm, allows us to reduce the scale from n x k down to poly log (n) x k. Once 
we reach the polylog(n) x k scale, we use a refinement from |2] of the classical canonical paths 
argument [12] . This allows one to obtain an upper bound on the relaxation time of a Markov chain 
in terms of the congestion ratio restricted to a subspace O' and the time the chain needs to visit O' 
with large probability. Here we use a further crucial input from |13j permitting us to identify a 
“canonical” subset of triangulations O' such that after T = 0{n^) the chain enters O' with large 
probability and such that the chain restricted to O' has small congestion ratio. A detailed high-level 


overview of the proof will be given in Section 4.1 


The rest of the paper is organized as follows. In Section we first recall some important tools 
from |3] and then formulate the main ingredients we need from [13]. Then, in Section]^ we develop 
the applications of improved canonical path techniques to our setting. In Section [^ we discuss the 
recursive scheme for the log-Sobolev inequality and prove Theorem o 


2 Main tools 

2.1 Triangulations with boundary conditions 

We will often consider subsets of n(n, k) consisting of triangulations in which some edges are kept 
fixed, or “frozen”; we call these constraint edges. Formally, let A' C A„^fc denote a subset of the 
midpoints, and fix a collection of non-crossing edges {r^, y G A'}, i.e., straight lines with midpoints 
in A' each of which connects two points of A[( ^ and passes through no other point of A]( ^. If 
a G Q{n,k) satisfies {ay = Ty , y £ A'}, we say that a is compatible with the constraint edges r. 
We interpret the constraint edges r as a boundary condition. 

We shall actually need a more general notion of boundary condition, in order to deal with the 
possibility of constraint edges whose midpoints lie outside the rectangle A° Let N be an integer 


3 



and consider the set Q%nk ~ • ■ • ,n + N} x {0,..., k}, i.e., a {2N+n) x k rectangle containing 

and let QN,n,k denote the set of midpoints of a triangulation of Q%nk- Fix a triangulation 
r of the region ^ ^ and call r the set of edges obtained from r by deleting some or all edges % 
with midpoint x G fc. Thus, r is a set of constraint edges for triangulations of such that 

all edges with midpoints in QN,n,k \^n,k are assigned. Given constraint edges r as above, we define 
n'^(n, k) as the set of all triangulations a of Q%nk fFat are compatible with r. Since the parameter 
N will play no essential role in what follows we often omit it from our notation. Since all elements 
of fe) have the same edges at midpoints in QN,n,k \ ^n,k, one can also view a triangulation 

a G Q'^{n,k) as an assignment of edges to midpoints in A^ ^ with certain constraints. Note that 
while the midpoint of a non-constraint edge of a triangulation a G k) is always contained in 

A„ fc, its endpoints need not be contained in A[^ we refer to Lemma 3.4 below for a quantitative 
statement on the smallest rectangle containing all non-constraint edges of any a G k) in terms 

of the length of the largest edge in r. 


The random triangulation a with boundary condition r is the random variable a G QJ{n,k) with 
distribution 

(3) 

where Z = k) We sometimes write n instead of /x’’ and instead of ^’^(n, k) if there 

is no need to stress the dependence on the constraint edges. We say that there is no boundary 
condition when A^ = 0 and the set of constraint edges r is empty. In this case Q^{n,k) coincides 
with n(n, k), the set of all triangulations of A° 


2.2 Ground states 

It is a fact that for any set of constraint edges r, the set of triangulations QZ{n,k) that are com¬ 
patible with r is non-empty. Among the compatible triangulations, we are particularly interested 
in those with minimal £i-edge length, which we call ground state triangulations. These are the tri¬ 
angulations of maximum weight in Q when A < 1, and they play a central role in our analysis. In 
the absence of boundary conditions, the ground state triangulations are trivial; every edge is either 
horizontal or vertical or a unit diagonal, so in particular the ground state is unique up to flipping 
of the unit diagonals. The presence of constraint edges can change the ground state considerably. 
However, the following result from [3l Lemma 3.4] reveals the strikingly simple structure of ground 
states for any set of contraints. 

Lemma 2.1. [Ground State Lemma] Given any set of constraint edges, the ground state triangu¬ 
lation is unique (up to possible flipping of unit diagonals), and can be constructed by placing each 
edge in its minimal length configuration consistent with the constraints, independent of the other 
edges. 

Given a set of constraint edges, we denote by d the unique ground state triangulation. (An ar¬ 
bitrary choice of the available unit diagonal orientations is understood in this notation.) If no 
confusion arises, we omit to specify the dependence on the constraint edges. An important struc¬ 
tural property of triangulations with constraint edges, which follows from Lemma |2.1[ is that from 
any triangulation a compatible with r one can reach the ground state a with a path in the flip 
graph with the property that no flip increases the length of an edge. 
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2.3 The Glauber dynamics 


The Glauber dynamics in the presence of a boundary condition r is defined as before (see equa¬ 
tion 0 ). with the modification that the midpoint x to be updated is picked uniformly at random 
among all midpoints of non-constraint edges. For any A > 0, this defines an irreducible Markov 
chain on k) that is reversible w.r.t. the stationary distribution fjT' (see [3] for details). It was 

shown in |3l Theorem 5.1] that for some constants C > 0 and Ai G (0,1), the mixing time of this 
chain in an n x /c rectangle satisfies Tmix < Ckn{k + n) uniformly in the choice of the constraint 
edges, whenever A < Ai. We also conjectured in [3] that the 0{kn{k + n)) mixing time should hold 
for all A G (0,1). 


2.4 Key ingredients from |13j 


We gather in Lemmas 2.2-2.5 below some estimates from m that will be crucial in our analysis; for 
the proofs see m Note that these estimates are valid throughout the sub-critical regime A G (0,1). 


The first lemma applies to the case where there are no constraint edges, so that the ground state is 
trivial. It follows from m Corollary 7.4], and establishes that after running the Markov chain for 
0{'n?‘) steps, the £i-length of a given edge has an exponential tail. For a given initial triangulation 
a = we denote by a* the triangulation after t steps of the chain. 


Lemma 2.2. Fix A G (0,1). There exist positive eonstants ci = ci(A) and C 2 = C2(A) such that 
for n > k > 1, for any t > cin^, any £ > 0, any midpoint x G An^k, o.nd any initial triangulation 
a G n(n, k): 

T (|o'^| > f) < Cl exp {—C 2 t)- 


The next lemma deals with the evolution in the presence of constraint edges r, and follows from m 
Theorem 7.3]. We denote by ax the ground state edge at x (compatible with r). Given a G n'^(n, k) 
and y G fc, we write Uy n Ua: / 0 if the edge ay crosses ax (not including the case where ay and 
ax intersect only at their endpoints). 

Lemma 2.3. Fix A G (0,1). There exist positive constants ci = ci(A) and C 2 = C 2 (A) such that 
the following holds for n > k > 1, for any set of constraint edges r. Let M be the ii length of the 
largest edge in any triangulation a G OF{n, k). Then, for any t > cikn{M + logn), and any £ > 0, 
we have 

IpfU n Ua, / 0} n > |aa,| < Ciexp (-C2^) . (4) 

Next we give a rough upper bound on the number of small edges intersecting a given ground state 
edge. We assume that a set of constraint edges r is given. For any triangulation a G FT{n, k), any 
ground state edge g, and any i G Z+, define 

Ig{a,i) = {cJa;, X G An^k- axFg and \ax\ < \g\ + £} . 

We denote by \Ig{a,tj\ the cardinality of Ig{a,i). For a proof of the lemma below, see m Propo¬ 
sition 4.4]. 

Lemma 2.4. Let g be a ground state edge, and let a G 12'^ (n, k) be a triangulation. 


i) If <^xF\ g ^ % then \ax\ > \g\, with strict inequality when the midpoint of g is not x. 
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ii) For any £ > 1, all midpoints of edges in Ig{a,i) are contained in the ball of radius 2i centered 
at the midpoint of g. 

Hi) There exists a universal c > 0 such that for any i > I we have 

\Ig(a,i)\ < ci'^ , and 1 < c£^. 

Finally, the lemma below establishes the probability of having a top-to-bottom crossing of unit 
verticals in a random triangulation a. By a “top-to-bottom crossing of unit verticals in cr” we 
mean a straight line of length k made up of k vertical edges in a each of length 1. The lemma 
below follows from m Theorems 8.1 and 8.2], 

Lemma 2.5. Let k G N and A G (0,1) be fixed. There exist positive constants c = c(A, k), 
6 = 6{X,k) and tuq = mo{\,k) such that the following holds. Let R be an m x k rectangle inside 
K,k with m > mo- Consider an arbitrary set of constraint edges r such that no edge from r 
intersects R. For any triangulation a G k), let Cji{a) be the number of disjoint top-to-bottom 

crossings of unit verticals from a that are inside R. Then, 

^{CR{a) <6m)< e-^^. 

Furthermore, let a, a' be two triangulations sampled from the stationary distribution p, given two 
different sets of constraint edges r, t' such that no edge of r, r' intersects R. Then, there exists 
a coupling of a, a' such that the probability that they have less than 5 m common top-to-bottom 
crossings of unit verticals is at most e~^^. 


3 Estimates via canonical paths 

We recall that the relaxation time Tj-ei is defined as the inverse of the spectral gap of the Markov 
chain. We start by showing that a direct application of the usual canonical path argument |12) 
yields an exponential bound on the relaxation time of the Markov chain that is valid for all A < 1. 
We recall the well known estimate relating Trei and T^i^ (see, e.g., [U Theorem 12.3]): 

Tmix < 7"rel(2 + log{l/p^)), (5) 


where = mina p{cr). 

Theorem 3.1. There exists a positive constant C such that for any X < 1, n,k G N and any set 
of constraint edges t, the Glauber dynamics on FF{n,k) satisfies 

Trei < exp(C'A:n). 

Before proving the above theorem we recall a useful structural fact. Given a set of constraint edges 
r and a midpoint x, consider the set Xlj. of possible values of ax, as a ranges in QF{n, k). Two edges 
ax,a'x G Flf are said to be neighbors if ax is flippable to a'x within some triangulation a G LF{n, k). 
Then it is known (see, e.g., m) that the induced graph with vertex set Llj. is a tree Qf. We will 
make use of the following technical lemma; see [3l Proposition 3.8] for the proof. 

Lemma 3.2. Fix a set of constraint edges r. For any midpoint x and any two triangulations 
a, a' G LF{n, k), the distance between a and a' in the flip graph is equal to Wx); where 

K.{ax,a'x) is the distance between ax and a'x in the tree Qf. 
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Proof of Theorem 13.11 For each pair a, a' G k), let be a shortest path between cr and 

a' in the flip graph. From Lemma 3.2, we have that for any triangulation r] in the path and 
any midpoint x, 

\'rix\ <\crx\y \(y'x\- ( 6 ) 


We can also assume that is a monotone path in the sense that it is composed of a sequence 
of edge-decreasing flips followed by a sequence of edge-increasing flips. 


Now, for any function /: —)■ M, we have 


a,al 


where we employ the notation = fiv) ~ fi'n')- For simplicity, below we write /x instead of 

and n instead of FF{n,k). Thus, using Cauchy-Schwarz, the variance of / with respect to /x 
satisfies 

Var(/) = \Y /^(^)/^(^')(/(o-) - 

77,77': 77^77' 


where p{r], r]') is the probability that the Glauber chain goes from r/ to rj' in one step, r] ^ rj' denotes 
that r] and rj' are adjacent triangulations, and we use the notation 




max 

77,77': 


E 

a,a': 


ip 

p{r])p{'n,r]'y 


( 8 ) 


for the so-called “congestion ratio.” Now assume that p{p,ri') > p{p',p), otherwise use reversibility 
to write p{ri)p{ri,ri') as p{ri')p{ri',rj). With this assumption we have that p{r],ri') > 2 |a^- Also, 
from Lemma 3.2 we have iFo-^o-'l = 0{nk{n + k)). The key property we use is that Q gives 


p{v) 

where we used the bound 


p { a ) p { a ') ^ ^-1 < ^-1 

X X 


Plugging this into Q, we obtain 

C(n) < Cnk{n + k) |n”(n, k)^. 


(9) 


Using Anclin’s bound [T] one has |U'”(n, A:)| < 2 A™.fcl. The proof is then concluded by recalling that 
Trei is the smallest constant 7 such that the inequality 

Var(/) < I Y Kr])piv,v')(yr!,r,'ff 

77,77': 


holds for all functions / : Fl'^ {n, k) 1 —)• M. 


□ 
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3.1 An improved canonical paths argument 

Here we establish a first polynomial bound on the relaxation time. The result here can be formulated 
as follows. 

Theorem 3.3. Fix A G (0,1) and k £ N. There exists a positive constant c = c(A, k) such that for 
any boundary condition r = {tx} such that \tx\ < n/4 for all x, the relaxation time of the Glauber 
ehain in k) satisfies 

Trel < n"". 


The strategy of the proof is as follows. We shah identify a subset of triangulations such that 
the congestion ratio C{Q') dehned as in but restricted to satishes a polynomial bound, in 
contrast with the exponential bound in ([9). Using a key input from m, we show that the Glauber 
chain enters the set H' with large probability after a burn-in time of T = 0(n‘^) steps. Following 
an idea already used in [2] we establish the desired upper bound on Trei by combining the above 
facts. 

We start with a deterministic estimate. 

Lemma 3.4. Let a G LF{n,k) he a triangulation of the n x k rectangle with boundary condition 
T = {tx} such that \tx\ < L for all x. Then, all edges of a are contained in the rectangle [—L,n + 
L] X [-L,k + L]. 

Proof. First, note that the ground state triangulation must satisfy the lemma, because all edges 
have size at most L. Now it is enough to show that there cannot be an increasing edge ax with 
X G such that a% [—L, n+L] x [—L, k+L] but ah edges of a are inside [—L, n+L] x [—L, k+L], 
We use the notation a^ to denote the triangulation obtained from a by flipping ax- In order to 
achieve a contradiction, assume that such an increasing edge ax exists and assume that a^ is at 
the left part of the triangulation (i.e., that its leftmost endpoint has horizontal coordinate smaller 
than —L). Let ay, az be the triangle containing ax such that the vertex v = ayCiaz has horizontal 
coordinate smaller than —L. Since a is completely inside [—L,n + L] x [—L,k + L], we obtain 
that ay and az are constraint edges. Also, since x £ A, ax must have one endpoint u of horizontal 
coordinate at least 0. This gives that ||u — u||i > L, and consequently, either ay or az has length 
larger than L, which is a contradiction. □ 


Next, we formulate a general upper bound on Trei in terms of the congestion ratio of a subset Ll' 
of the state space H, a time T, and the probability needed to reach Ll' within time T. A version of 
this lemma appears in [21 Theorem 2.4]. For the reader’s convenience we give a detailed proof. 

Lemma 3.5 (Canonical paths with burn-in time). Consider a Markov chain with state space 
irredueible transition matrix p{-, •) and reversible probability measure pL. Let Ll' C Ll be a subset so 
that between eaeh a, a' £ Ll' there is a path Fo-^g-/ in the Markov ehain that is entirely contained in 
Ll'. Define the congestion ratio 


c{n') 


max 


E 

cr,cr': 


^i{a)^l{a')\T^^„'\ 

h{ri)p{v,V') 


( 10 ) 


where the sum is over all pairs of states a, a' £ so that the path Ta,cF' uses the transition (r/, ??'). 
Fix T G N and let p he a lower bound on the probability that at time T the chain is inside Ll', 




uniformly over the starting state in O. Then the relaxation time satisfies 

^ 3Cin') 


Trel <-1" 


P P" 


Proof. We run the Markov chain for T steps. For a, r G O, let ^o-(r) be the probability that, 
starting from a, the Markov chain is at r after T steps. Note that /Uo-(^0 — P- cr,T G Q, and 
for any path 7 of length T in the chain starting at a and ending at r, let be the conditional 

probability that, given the initial state a at time 0 and the final state r after T steps, the Markov 
chain traverses the path 7 . Then, for any function / : 17 —)• M, we have 

Var(/) = | E = \ E E ^ 


cr,cr'Gf^ 


cr,cr'Gr^ r],7]'£Q' 


X + See72 + Seer„_„, Ve/) , 

71)72 

where the three sums inside the parenthesis are over the edges of the paths 71 , 72 , and Tri,ri'-, 
respectively. Then, applying Cauchy-Schwarz, we obtain 

Var(/)<2 E E b,(n')M..(n-) 

cr,(T'eo»?,7e!7' ^ ’ 

X Z^.,,(7l)^.',7(72)(rEee7l(^‘=/)' + ^Eee72(^e/)2 + |r,,7| Eeer,.,,(Ve/)2). 


We write the right-hand side above as yli + yl 2 + , where 




(T,cr'Gr2 T],r]'£Q^ 


/N Pa{v)Pa’{v' 


^ W. .VO ^co,o 


(T,cr'Gr2 Tj^rj'GQ' 


^3 = 2 Z] 


o-,o-'en ri,r]'&Q' 


^ etr2 


We start with vli. Summing over 72 ,( 7 ', 7 ', and using na-'^n'A) = 1, we have 

o-enrjeO' '71 ee 7 i 

Changing the order of the summations, and summing first over all pairs of adjacent states r ~ r'. 
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we get 


Al = -T ^J-iT)p{T,T'){Vr,T'ff Y 


n{a)pa{ri)va,ri{l) 


< 


< 


‘T-p 

^ KT)p{'r,r'){Vryf)" 

: T'^t' 


3T 


po po/ r 
a£U,ri£U','y: (r,T')E 7 

p{a)pa{p)l^a,'q{l) 


E 


(TGf^,77Gr^',7: (t5t')G7 


p{t)p{t, r') 


< 


2p 

3T 


2p 


E , . , .. 2 1^44 (Ms-rkov chain traverses (r, 7) within T steps) 

p{T)p{T,T){Vryf) --7T- 

^(r)p(r,rO 

> ^(r p(r,r (v^y/ —— -— = - V{f,f), 

^ p{r)p{T,r) p 


where denotes the measure induced by the Markov chain started from stationarity, and we 

use the notation 

^(/>/) = ^ Y P{T)p{Ty'){Vryff (11) 

T,r'Gr2: r~r' 

for the so-called Dirichlet form. For the second term, we have by symmetry that A 2 = Ai. For ^ 3 , 
we use p < and sum over 71 , 72 , cr, a' to obtain 

E P{ri)p{ri')\T^y\ Y 

riy&w 


Changing the order of summations, we get 


As<^Y^V{fJ). 


The result now follows since Tj-ei is the smallest constant 7 such that the inequality 


Var(/) <7iP(/,/) 


holds for all functions / : fl 1 —)• 


□ 


Proof of Theorem 13.31 Let T = cin^k for some large enough constant ci = ci(A) > 0. Thanks 
to Lemma 3.4 we may apply Lemma 2.3 with M = 2n + k. Thus, for any given x G A„ ^ and 
ground-state edge ax with midpoint x, taking £ = C 2 log|A„^fc| for some large enough constant 
C 2 = C 2 (A) > 0, and taking the union bound over all x G A„ ^ in Q we obtain that the triangulation 


cr 


at time T, for an arbitrary initial condition a, satisfies 


(T n <4. 0} n {ly > |5,| + <}) < 


-1 


n 


( 12 ) 


Let 


7 = |cr: for all x,y £ Any, |cra,| < Icrj,! -|- i and 1 {ay n cJa, / 0 ) < 1 (|(Tjy| < |cr 3 ;| -|- ^) |. 


Thus (12) implies that P (u^ G 7) > 1 — n Note that fl' is a decreasing set in the sense that if 
cr G 7 then for all cr' that can be obtained from a by performing decreasing flips, we have cr' G 11'. 
This allows us to construct a path F^y within 11' between any pair of triangulations cr, a' G H'. 
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We now describe the path Fix two triangulations a, a' G O', and any midpoint x G 

Let g' be a ground state edge at x. The edges that need to be flipped to transform ax into a'^. are 

we have 


contained in Ig{a, i) U Ig{(x', i) (recall the definition of Ig from Lemma 2.4). By Lemma 


2.4 


that all edges in Uo-en' have midpoint inside a ball of radius 2i centered at x. This implies 

that if we partition [0, n] x [0, k] into slabs of horizontal width 2i, we can find a sequence of flips 
that transform a into u' slab by slab, from left to right, so that when transforming the zth slab, 
only edges with midpoints in the ith and {i + l)th slabs need to be flipped. In each slab, we just 
perform the minimum number of flips needed to transform that slab into cr', and we do that by 
first performing all decreasing flips and then all increasing flips. 


Our goal is to apply Lemma 3.5, for which we need to bound the value of the congestion ratio 
C(n'). To do this, consider a pair of adjacent triangulations r],r]'. Assume that r/, ry' differ at an 
edge of the fth slab. Therefore, if a, a' are two triangulations for which the path between them 
includes the transition (rj, ry') we know that triangulation ry has slabs 1, 2 ,..., i — 2 equal to a' and 
slabs i + 2,i + 3,... equal to a. Let ^ be a partial triangulation in Q' of the first i — 2 slabs and 
m be a partial triangulation in fl' of the middle slabs so that m and a are compatible, meaning 
that m and the edges of a inside slabs i + 2,i + 3,... can coexist to form a full triangulation. 
Similarly, let be a partial triangulation in 12' of the last slabs (z + 2,i + 3,...) and m! be a partial 
triangulation in 12' of the middle slabs so that ^', m! and the edges of a' inside slabs 1, 2 ,..., i — 2 
are compatible. Assume that p{r],r]') > p(rj',rj), which implies that p{T],r]') > (otherwise, 

replace p{r])p{r],'n') with p{ri')p{ri'iP) in C(f2')). Let pi be the part of ry inside slabs r — l,i,r + 1. 
Then, summing over all ^', m, m' as above such that (ry, ry') is a transition in the path from m, a 
to a\m', and noting that the path between a and a' has length at most 22'|A„fc|, we obtain the 
following upper bound for C(I2'): 


c{n')<u\K,k? Y. 


Zn' 


where Z^i = Ylaen' Instead of summing over m, m', we will sum over triangulations m" of the 
middle slabs that are compatible with both ^ and and are to be interpreted as m A m'. Given 
m", we sum over rrr, m! that can be obtained from m" by increasing flips and such that (ry, ry') is a 
transition in the path from m, a to ex', rrr', ^'. Let A(m", rrr, rrr', ry) be the indicator that all four of 
them are compatible, as described above. When A{m", rrr, rrr', ry) = 1 we have that |rya;| < |rrr 3 ;| V|rrr(,| 
for any midpoint x in the middle slabs. Hence, |rrr"| + |rrr \ rrr"| + |rrr' \ rrr"| > |ryj|, which gives 


C{Q') < Ai\An,k\ 


' E 


Aki+k'i+i"*"i 


E 

m,m‘ 


+ ry). 


Since A < 1, we can simply use Anclin’s bound [T] saying that the number of triangulations of an 
I X k region with arbitrary constraint edges is at most 2^^^ to obtain that 


C(f2') < 4^|A„,fcr2 


2 r) 6 /c^ ^ ^ 


Aki+k'i+|m"i 


< U\K,k?2^'^^. 


Plugging everything into Lemma 3.5 completes the proof. 


□ 
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4 Proof of Theorem 11.11 

4.1 High-level overview 


The proof is composed of three main ingredients: (i) a good ensemble, (ii) a decay of correlation 
analysis, and (iii) a recursion for the logarithmic Sobolev inequality. 


The good ensemble. The first step is to show that uniformly over the initial condition, with high 
probability, for all times t £ [T,T + n^], with T = 0{n^), the Markov chain stays within a subset 
of triangulations where all edges have length at most Clogn for some constant C > 0. We 
will call this subset the good ensemble. This result will be a consequence of the tail estimate of 
Lemma 2.2 Therefore, we will couple our evolution in the time interval t £ [T, T + n^] with the 
Markov chain restricted to the good ensemble, which evolves as before, by attempting to flip edges 
chosen uniformly at random, but with the suppression of any edge flip that would render an edge 
longer than Clogn. The structural properties of triangulations imply that this Markov chain is 
irreducible. Moreover, the reversible probability measure is given by ^ = /u(- | 0), the measure n 
conditioned on the event a £0.. Since /r and Ji can be coupled with high probability, it is sufficient 
to analyze convergence to equilibrium for the restricted chain, and to show that the latter mixes in 
time T' = 0{n^). We will actually prove that the restricted chain mixes in time T' = npolylog(n). 
For the rest of this discussion we assume that we are working with the Markov chain restricted to 
the good ensemble 12. 


Decay of correlations. We split the set of midpoints ^ into two intersecting slabs and Ar, 
where A^ contains all midpoints with horizontal coordinate smaller than n/2 + 2(7 log n and Aj- 
contains all midpoints with horizontal coordinate at least n/2 — 2(7 log n. Note that A^ n Aj- is 
a slab of height k and horizontal width 4Clogn. Let be the cr-algebras generated by the 

edges with midpoints in Aj. \ A^, A^ \ A^. respectively. We want to show that, conditional on any 
event F £ the distribution of the edges in A^ \ Ar is not affected much, and similarly for events 
F £ Ti. The intuition for this is that the intersection A^ n Ar of the slabs is large enough to allow 
correlations from A^ \ Ar to decay. We will make this intuition rigorous by showing that there exists 
a positive e = e(A) such that, for all J^-measurable functions fi and all TT-measurable functions 
/r, we have 


sup |^(/£ I F) -/i(/£)| < n ^ll/^lli and sup |;u(/r | F) - //(/r)| < n ^||/r||i, (13) 


where /x(/ 
Lf 


to denote the 


norm 


F) stands for the expectation of / given the event F and we use 
1 = 

The high-level argument for ( |13[ ) is the following. Fix any valid collection of edges with midpoints 
in A^ \ Ar, that is, a partial triangulation from fl. This defines an event F £ Fi. We will construct 
a coupling of one triangulation a distributed according to ]!{■ \ F) and another triangulation a' 
distributed according to /!(•). We do this by first sampling the edges of a' whose midpoint is 
in A^ \ Ar. Call this event F' £ Fi. Since we are restricted to the good ensemble, the edges 
of F and F' have length at most (71ogn. Therefore, none of them crosses into the right half of 
A^n Ar. Lemma 2.5 therefore ensures that we may couple the sampling of edges in Ar so that, with 
probability at least 1 — we put the same top-to-bottom crossing of unit verticals in a and 

a' inside the right half of A^ (7 Ar. In particular, this implies that we can couple a and a' so that 
they agree on Ar \ A^. This will establish (13). 


The log-Sobolev inequality. An important ingredient in the proof of Theorem o is the use of the 
logarithmic Sobolev inequality for the good ensemble. For any positive function /, let Jl{f) stand 
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for the expectation of / in the good ensemble, and let 

Ent(/) = ^ (/ log log 

cr 

denote the entropy of /. Also, define 

(T,(T^Gr2 


where 


p{a, a') = 


Ak'l 


AkI + 


1 (cr ~ fj') 


As usual a ^ a' means that a, a' differ by a single edge flip. Note that p{a,a') = |A„^fc|p((T, cj'), 
where p is the transition matrix of the discrete time chain. Thus £{f,f) can be interpreted as 
the Dirichlet form of the continuous time Markov chain where every edge of the triangulation 
independently attempts to flip at rate 1. 


Let cs be the log-Sobolev constant of this Markov chain, defined as the smallest constant c > 0 such 
that for all functions / one has 

Ent{f)<c£{f,f). (14) 

It is known (see e.g. m Theorem 2.9]) that cs is related to the mixing and relaxation times via 


^ mix 


< ^ (4 + log+ log and 2 T^ei < cs < T^ei ^ 


(15) 


where = min^^^ ^(cr), and we use TrepTmix to denote the relaxation time and the mixing time 
of the continuous time chain restricted to the good set. These bounds should be compared with 
Q. In particular, it will be crucial for us to work with the log-Sobolev constant rather than the 


relaxation time in order to obtain the strong bound on mixing time claimed in Theorem 1.1 


Recursion. We will bound the (restricted) log-Sobolev constant via the so-called bisection method 


introduced in m- Let A£, Ar and be as above. Using the decay of correlations in the 

decomposition estimate in 0 Proposition 2.1] implies that for all functions / : II i—)> M we have 


Ent{f)<{l + 0{n-^))p[Ent{f\Te) + Ent{f\T,)] < {l + 0{n-^)) 2c^^^£{fJ), (16) 

where Cg^^ is the largest log-Sobolev constant among the systems conditioned on J-i and J-r and 
the factor 2 comes from the double counting of flips within the region A^ n A^. Hence, we obtain 
that Cs < (1 -|- 0{n~^)) 2cg^\ We would then like to recursively apply the same strategy to bound 
Ent(/^ I T^) and Ent(/^ | T))- Indeed, p{- \ J-)) is a Gibbs measure on triangulations with 
midpoints in A^, and we may split Ai into two intersecting slabs, establish decay of correlations 
and again use the decomposition above to further reduce the original scale. One caveat is that now 
we have to take into account the boundary conditions dictated by the conditioning on T). These 
consist of constraint edges protruding from the right boundary, with midpoints in Aj. \ A^. The 
boundary conditions will not be a major problem since we are in the good ensemble so these edges 
cannot protrude more than a distance C log n. After j such iterations, we will be considering slabs 
of size roughly n2~^, with edges of size at most Clogn protruding from both the left and right 
boundaries. It will be convenient to iterate this procedure for j = j* steps, where n2~^* is roughly 
log® n, so that protruding boundary edges are still far away from the middle of the slab, which is 
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the crucial region for exploiting the decay of correlations. With this strategy, after j* iterations we 
obtain 


cs < {i+o{n-^)y* 


Employing the general polynomial bound on the relaxation time of Theorem 3.3 and the relation 
between cs and T^ei, we obtain that c^*'^ is at most polylog(n) uniformly over all boundary conditions 
in the good ensemble. The main problem is that the term 2^* is too large (of order ^ by our 
choice of ). As in |10j we overcome this difficulty by randomizing the location of the split of A„ ^ 
into Ai and Aj-, and similarly for the other scales. The idea is to first split into three disjoint 
slabs with height k, the left and right slabs with horizontal length |(n — log^n), and the middle 
slab with horizontal length log^ n. Then we further split the middle slab into smaller slabs (that 
we call rectangles) each with horizontal length 4(7 log n. We choose one such rectangle uniformly 
at random, and define A^ to be the midpoints to the left of this rectangle (including the rectangle) 
and Ar to be the midpoints to the right of this rectangle (including the rectangle). With this 
randomization, (16) will be improved to 


Ent(/^) < {l + 0{l/log^n))cy'^S{f,f), 


where log^ n is roughly the number of rectangles in the middle slab of An^k- Then, iterating j* 
times (with j* as above) we get 


Cs < (1 + 0(j*/log^ n)) = polylog(n). (17) 

Once we obtain , using ( |15[) we can conclude that the continuous time Markov chain restricted 
to the good ensemble satisfies Tmw = polylog(n). From this the desired conclusion for the discrete 
time Glauber dynamics will follow in a simple way. 

We now proceed with the detailed proof of Theorem EH 


4.2 The good ensemble 

Let (T®, cr^,... be the discrete time Markov chain on triangulations of A® ^ with no constraint edges. 
The first step is to show that after a burn-in time of order n^, during a very long time interval, 
the largest edge of the triangulation is of order at most logn. Let C = (7(A) be a large enough 
constant, and dehne 

O = |(T G 0: {o'xl < (7 log re for all x G (18) 

The set represents the good ensemble. The next lemma will allow us to analyze the Markov 
chain restricted to the set 11. 

Lemma 4.1. Fix A G (0,1). There exists a constant ci = ci(A) so that if we set T = cire^ then for 
all n> k > 1 

•* (nK/ w‘ £«})>!- 


Proof. For any given x G A^^k and any t > cire^, Lemma 2.2 gives that 


P (|(T* I > (7 log re) < exp(—C 2 ( 71 ogre), 


for some constant C 2 independent of (7 and re. Setting (7 large enough and taking a union bound 
over all x G An^k and all integers t G [T, T -|- re^] concludes the proof. □ 
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4.3 Decay of correlations 

Let r C A be a slab of width w; that is, for some x G Z, 


r = An,k r][x,x + w] X [0,k]. 

We assume throughout that w > ^ C^log^n, where C is fixed as in (18). 

Partition P into three slabs, two of width roughly — C^log^n) and one of width roughly 
log3 n. More precisely, for P as above, let 

= + x[o,fc], r2 = A„,fcn(x + ^^^=^^,x + ^^^±^^] x[o,fc] 

and Ps = An^k H (x + ^ + rc] x [0, A:]. 

Partition the middle slab r 2 into disjoint slabs Ji, J 2 , ■ ■ ■, Jg (from left to right) each of width 
dClogn, with 


s = 


log3 n log 2 


n 


4C log n 4 

Let i be an integer chosen uniformly at random from |l,2,...,s}. Finally, define 
P^ = Pi U Ji U J 2 U • • • U J, and Pr = P 3 U U J,+i U • • • U 


(19) 

( 20 ) 


Then, P^ represents the left portion of P, Pr represents the right portion of P, and P^ n Pr = J^. 

We need to introduce some more notation to be precise about boundary conditions. For any u G fl, 
A C A„^fc, if (T = {ax, X G An^k} then we write a a for the set of edges {ax , x G A}. If ^ = ua for 
some (T G n and A C A^^k we say that a contains ^ and we call ^ a partial triangulation in 11. If 
Ar\ A' = % and ^ = a a, = cta’ for some u G H, then we define ^ U ^' = <tauA'- 


We use partial triangulations in 11 as boundary conditions for a region B <Z T. Fix a partial 
triangulation We denote by A^ C An^k the set of midpoints of the edges in Let 11^ denote 
the set of full triangulations cr G H that contain We define for any C P, and any ^ such that 
Ag C An,k \B, 

= {as '■ (T G n^}. (21) 


For any rjB ^ H^, let 


a|(^b) = 


E 


CTgOS: crg=riB 




/i(ll€) 


be the induced probability measure over ll|j. In words, /i|j is the marginal distribution over 
midpoints B when we impose a boundary condition If ^ is empty (no boundary condition) we 
simply write and ^b- 


Lemma 4.2. There exists a positive constant c = c(A, k) such that for any partial triangulation 
f with Ag C An^k \ L, for all functions /^,/r : H 1—M such that fi depends only on edges with 
midpoint m P^ \ and f„ depends only on edges with midpoint in Pr \ and for any a^ G 
and fJr G Up y , we have 

and 
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Proof. We will establish only the hrst estimate; the second follows by a symmetrical argument. 
Since fi depends only on edges with midpoint in r£\ it is enough to show that, for any cJr G rip y 

and any r G , we have 



< exp(-C2Clogn), 


( 22 ) 


for some positive C2 = C2(A, k), where C is the constant in the dehnition of the width of 


Let r] and rj' be random triangulations distributed as and respectively. Let P denote 

the following coupling between rj and r/'; refer to Figure The idea is to sample recursively edges 
from the pair ( 77 , rj') in vertical strips inside from right to left from a suitable coupling of //j and 

to ensure that, with large probability, there is 


k-t 


Here we will use the estimate of Lemma 


2.5 


a common top-to-bottom crossing of unit verticals within J^. On this event we can safely resample 
iVre\j,,riP\j) in such a way that = r. 


We now present the details. Consider the midpoints of F in order of their horizontal coordinate, 
from largest to smallest (i.e., from right to left in Figure]^. Let vq be the leftmost integer horizontal 
coordinate of points in Fr \ J^, and let Vq = ? U Ur and Hq = ^. Now for i > 1, dehne Uj, V), V)' 
inductively as follows. Let Vi < Vi-i be the rightmost integer horizontal coordinate that is not 
crossed by an edge of Vq U Vi U U V 2 U H 2 U • • • U Vi-i U Vf_i. Using the coupling from Lemma 
sample all edges of r] and rj' whose midpoints have horizontal coordinate Uj, and denote them by Vi 
and Vf, respectively. There are two cases. In the hrst case, at least one edge of V) or Vf is not a 
unit vertical (as happens with i = 1, 2 and 3 in Figure]^. In this case, continue by dehning Uj+i as 
described above. If Uj+i is a horizontal coordinate in J^, sample U+i and VPi as described above 
and iterate. Otherwise, if Uj+i is not in J^, stop this procedure and sample the remaining edges of rj 
and rj' independently. In the second case, all edges in V) and Vf are unit verticals (i.e., they create 
a top-to-bottom crossing of F, as in Figurefor 7 = 4). Then stop the procedure above and sample 
the edges with horizontal coordinate smaller than Vi identically in both 77 and 77 ' (as depicted by 
the gray edges in Figure]^, and then sample the remaining edges (that necessarily have midpoints 
in Fr) independently in 77 and 77 '. Let be the event that 77 and rj' have a common top-to-bottom 
crossing of unit verticals with midpoint in J^. 


2.5 


Let 77 ^, 77 ^ be the edges of 77 , 77 ' with midpoints in F^ \ J^, and let rj^,ri[ be the edges of 77 , 77 ' with 
midpoints in Fr \ J^. Using the above coupling, for any r' G ^p\j we obtain 


= r,r]g = r') 

r^\Jt 

= P( 77 ^ = r', 77 ; = r') ^ F{r]e = T,r]'^ = r'). 

The hrst term on the right-hand side above is at most IP( 77 ^ = r') = Mr)'\j second term is 

bounded above by 


^(77^ = T')IP(r 7 (; = r') < = 


where the last step follows from Lemma 2.5 
the terms, we obtain 

AUCT; 




I = t') < = r')exp(- 4 cC'logT 7 ), 

Plugging this into the equation above, and rearranging 
exp(- 4 cC'log 77 )) 
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re\J, 


rr\ J. 


Vi 


V3 V2 


Vl 


Vo 




Figure 2: Coupling between (above) and (below). Note that the figure is not to scale: i 

reality, the middle region is much smaller than the two outer regions. 


m 


which holds uniformly over r' and Similarly, we write 

^ = r)+ P(% = \ m = r) 

- ^rAJ.exp(-4cC'logn). 


and the proof of (22) is completed by rearranging the terms and setting C 2 appropriately. 


□ 


4.4 Recursion via bisection 


We consider slabs of different scales: we index the scale by j, where j = 0 corresponds to the full 
slab of width n, while at scale j, we have slabs of width w roughly equal to n2~K The finest 
scale will be 

j* = min {j > 0: n2~^ < C® log® n}; (23) 

in particular, n2~^* > ^(C^log^n). Recall how slabs are split and the definition of i from the 


construction of F^ and Fr in the paragraph culminating in (20). 


Consider a given scale j G {0,..., j*}, and let F = Fj be a slab at scale j. Set Wq = n, and define 
the intervals 

Wj = [n2~^ - jC^ log^ n, n2~^ + jC^ log^ n] , j = 1,..., j* 

Notice that our slab F is obtained after j steps of the bisection procedure, so that F necessarily has 
width w G Wj. Let cr G 12 be an arbitrary triangulation in the good ensemble and set ^ ^ 

12fe\r a- boundary condition for the region F. Consider the continuous time Markov chain on 


12p with Dirichlet form 


4(/./) = 2 


OT,(TpGnp 


Ar('^r)/5r(oT, (Tr)(/(crp U ^ - /(c^r U 0)^ 


(24) 


where / : 12 i—)• M and 




( 25 ) 
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Let cs(r,(^) denote the log-Sobolev constant defined as the smallest constant c > 0 such that 

En4(/2)<c4(/,/), (26) 

holds for all functions /, where Entp(/^) denotes the entropy of with respect to /Xp. 

Einally we define, for each j, 

7 j = sup {cs(r, ^): r C An^k is a slab of width w G Wj, and ^ G fc\r}' 

The following lemma summarizes the result of this recursion. 

Lemma 4.3. There exists a positive constant C 2 such that, for any integer j G {0,... ,j* — 1}, 

7, <(l (l + 7,.,,. 


Proof. Let E be a fixed slab of width w G Wj, and let ^ be a given boundary condition. Let s, i, L^ 
and Lr be as described in the paragraph culminating in (20). From Lemma 4.2 and [5l Proposition 
2.1], for any function /: 12 i—>■ M we have that Entp(/^) is bounded above by 

( 




t=i 


„-C 2 logn 


(27) 


?U(T^ / f2\ 






Note that Entp'^'^’’ (/^) and Entp^*^^ (/^) are entropy functions for slabs on scale j + 1 given boundary 
conditions ^ U Ur and fUae, respectively. Therefore, by (1^ we have 


< cs(r,,^UcT7)47^(/,/) < 


?Uo-r , 


7?Ucrr , 


(28) 


and similarly for the second term in (27). Now we claim that 




i=i 








To prove (29) we proceed as follows. Since a given edge ax in a triangulation has at most one value 
(t(, / ax it can flip to, we may write the flip rates (25) as 


pf((Tr,cJr) = 


Al-il 


xer 


Ak^l + aKI 


(ar ~ 0-^ U^; ax ^ a!,) =: Y^^iri^r) ■ 

xer 


Therefore, 


^kf’f) = ^Y Y f^r(^r)pir(^r)(Vxf(arUCf, 


(30) 


orsn^ 


where we use Vxf to denote the difference in values of / before and after the flip at x. It follows 
that 


Yj ^rr\j,('^r)<f^r^ "^(/j/) 

^\Y Y Y 4r'(^rJpi77(^r,)(Vx/(cJ7,j,^ U^U(Tr))^ 

a:er« co^^'^r 
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where, as before, we use the shortcut notation Using 

= Ar(^U U f^r)4,r(^U U cTr) 

and rearranging the sum, we obtain 




o-psof 


A similar expression holds for the second term on the left-hand side of (29), and the desired estimate 


follows from the expression (30). 


Plugging (29) and (28) into the bound in (27) we have 


Ent' 


« if) < (l + e-'^^ 7i+i (1 + i) 4(/> /)• 


This establishes that cs(r,^) < (l -|- e '^aiogn^ Ij+i (l + ^)- Since this bound does not depend on 
^ and the choice of slab T at scale j, the proof is completed by using the value of s from (19). □ 


We conclude the proof with the base of the induction. 

Lemma 4.4. There exists a constant c = c(A, k) such that 

7A < log'^ n. 

Proof. Let T be a slab at scale j*, so that the width of P is of order log®n. Let ^ G ^Ar,k\r 
be a boundary condition. We note that the argument of Theorem 3.3 can be repeated with no 
modifications for the chain restricted to the good set U. Therefore, there exists a constant ci = 
Cl (A, k) independent of T and f such that the relaxation time of the discrete time chain on T with 
boundary condition f is at most log'^^ n. Passing to continuous time, we have that T].ei(r, ^ < 
log'^^ n. Since triangulations in Up have edges of length at most Clogn, there exists a constant C 2 
such that 

min Ar('^r) > n~^^, 
orsfip 

uniformly over all slabs P at scale and boundary conditions Therefore, using the relation 


between the relaxation time and the log-Sobolev constant from (14) we have that 


cs(r,6<7;ei(r,o 


log(n'^' 

1/2 


Since the bound above is uniform T and this proves the desired bound on 


□ 


4.5 Completing the proof 


Proof of Theorem 
on U. Lemma 4.3 


1.1. We start by bounding the mixing time of the discrete time Markov chain 


implies that the log-Sobolev constant of the continuous time Markov chain on 


An^k with no boundary condition is at most 


cslAnU < 70 < (l + (l + cA?;;)’'''< 27i 


vi*-i 


19 



















where the last step follows since j* < log 2 n. Also, we have that 


min /r(cT) > 

ctEH 


_\|A„,fe|Clogn 


where comes from Anclin’s bonnd of for the number of lattice triangulations [T], 

and the fact that the total edge length of any triangulation is at least |A„^fc|. Therefore, using the 


relation between the mixing time and log-Sobolev constant in (14), we deduce that the mixing time 
2hiix of the continuous time Markov chain on is bounded above by logn. Thus, the mixing 
time of the discrete chain in H is at most | A„ log n, for some constant c. Using Lemma 4.4 and 
the fact that |A„^fc| is of order nk, we obtain that the mixing time of the Markov chain restricted 
to n is at most cnlog'^n, for some new positive constant c (which depends on k and A). 

Now we compare the restricted chain on 11 to the original unrestricted chain on H = ll(re, fc). 
Let Ti = cnlog'^n and fix the constant c > 0 so that the total variation distance between the 
restricted chain at time Ti and the restricted stationary distribution J1 is at most 1/8. We obtain 
the mixing time of the unrestricted chain via the following coupling. Let Tq = cire^, where ci 
is the constant in Lemma 4.1 Let the unrestricted Markov chain run for Tq + Ti steps. With 
probability at least 1 — n“^, the unrestricted chain never leaves the set H during the time interval 
[ro,To + Ti]; therefore, we can couple its steps with those of the restricted chain. This gives that 
the total variation distance between the unrestricted chain at time Tq + Ti and the stationary 
distribution is at most + 1/8 + /u(H \ H). Sin ce H \ H only contains triangulations for which 


the largest edge is larger than Clogn, Lemma 2.2 ensures that /u(H \ H) < n ^ for large enough C, 


and therefore the total variation distance between the unrestricted chain at time Tq + Ti and its 
stationary distribution is at most 1/4. This completes the proof of Theorem □ 
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